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Persistent current in an almost staggered Harper model 
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In this paper we study the persistent current (PC) of a staggered Harper model, close to the half¬ 
filling. The Harper model, which is a quasi-periodic system, is different than other one dimensional 
systems with uncorrelated disorder in the fact that it can be in the metallic regime. Nevertheless, 
the PC for a wide range of parameters of the Harper model does not show typical metallic behavior, 
although the system is in the metallic regime. This is a result of the nature of the central band states, 
which are a hybridization of Gaussian states localized in superlattice points. When the superlattice 
is not commensurate with the system length, the PC behaves as in an insulator. Thus even in the 
metallic regime a typical finite Harper model may exhibit a PC expected from an insulator. 

PACS numbers: 73.23.Ra, 71.23.Ft, 73.21.Hb 


In the last decade, quasi-periodic one-dimensional (ID) 
potential, also known as the Aubry-Andre-Harper (for 
short Harper) model (THU, has garnered much interest. 
The main reason for this interest is that it is the dis¬ 
order of choice for cold atoms, since it may be created 
by the superposition of two incommensurate periodic po¬ 
tentials [5- m Contrary to white noise disordered sys¬ 
tems which are always localized for ID the Harper 
model shows a ID metal-insulator (Anderson) transition 
for non-interacting particles as function of the strength of 


the potential [2H2j]. This has been demonstrated experi¬ 
mentally in cold atoms , as well as for optical systems 
[§, 9j. Additional effort has gone into understanding the 
influence of electron-electron interactions on this metal- 
insulator transition 13 and on the the many-body 
localization in these quasiperiodic system T(| 17|. The 
Harper model also exhibits topological edge state [3, 0, 
and show counterintuitive behavior of the compressibility 

ME|. 


The current carried by the ground state of a ring 
threaded by a magnetic flux is called the persistent cur¬ 
rent (PC) [22|. In the presence of elastic scatterers (dis¬ 
order), the persistent current is suppressed H-Eljj. In 
the diffusive regime, on the average it is suppressed by 
a power law of the circumference of the ring L, while 
in the localized regime it suppressed exponentially by 
exp (—L/C), where £ is the localization length. Since for 
non correlated (white noise) disorder, all states of a ID 
system are localized, one expects that for ID systems the 
averaged PC is always suppressed exponential. 

The PC is a very effective way to evaluate the sen¬ 
sitivity of the system to boundary conditions, i.e., the 
conductance of the system ^3, ^ , and therefore a great 
way to identify whether your system is metallic or local¬ 
ized. Thus, naively, we would expect that the PC for the 
Harper model in the metallic regime (i.e., not too strong 
on-site potentials) will on the average be only weakly 
suppressed by the potential. Here, we will show that 
the persistent current of the Harper model can exhibit 
a rather intricate behavior, which can skew the simple 
picture presented above. 


In this paper, we study the PC of a Harper model of 
spinless fermions on a ring threaded by a magnetic flux. 
This is a tight-binding model in which the on-site poten¬ 
tial is spatially modulated with an irrational frequency. 
We would focus here on irrational frequencies which their 
modulus is close to half. This corresponds to a fast (two 
site) modulation with a slow envelope. These frequen¬ 
cies exhibit an increase of the compressibility when the 
electron-electron interactions are increased, opposite to 
the influence of interactions in regular disordered systems 
M l2ll | . Here we will show that for this range of frequen¬ 


cies close to half-filling, the PC shows a non-monotonous 
dependence on the systems size, where for most values of 
L the PC is strongly suppressed. 

The tight-binding Harper model Hamiltonian for spin¬ 
less fermions on a ring threaded by a magnetic flux is: 


l l- 1 

H = kcos(2nbj + </>)cjcj — t (e^^cjcj+i + /i(d) 
0 =1 3 =1 

where Cj is the single particle annihilation operator 
on site j, t is a real hopping amplitude. The magnetic 
flux is denoted by <p, and (po is the quantum flux quanta 
Po = hc/e. The strength of the on-site potential is con¬ 
trolled by A > 0. The on-site potential is modulated by a 
frequency b , and <fi is an arbitrary phase factor. It should 
be clear that since we are interested in a ring, <f> is irrel¬ 
evant and will be ignored through the rest of this paper. 
We will be interested in the metallic regime of the model, 
i.e A < 2 1. The irrational frequency may to written as 
b = Z + 1/2 + e, where e is irrational. Therefore, we can 
write the on-site potential as: 


cos(27r6j) = cos(27tZj + n j + 27rej) (2) 

= (— iy cos(27rej) 

When e 1/2 the system is called an almost stag¬ 
gered Harper model for which the fast frequency of the 
(—1) J term is modulated by the slow frequency, e, of the 
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cos(27rej) term. In the almost staggered case the energy 
spectrum of system shows unique features such as a cen¬ 
tral band that is separated from the other bands by a 
large gap, of order A, as can be seen in FigQ] Also the 
two bands sandwiching the central band show similar fea¬ 
tures, i.e., a rather narrow (flat) band and a large gap to 
the next band. Changing the length of the sample length 
from L = 900 to L = 1000 does not change its gross 
features, although some difference in the energies of the 
edg e states in the gaps are apparent. As detailed in Ref. 
[20], for e <C 1, there are L n = 2|e| L states in the cen¬ 
tral band, corresponding to the number of intersections 
with zero of the slow modulation envelope, which occur 
at cos(27 rejVi) = 0. These valleys are shown in Fig(2]for 
a smaller system (L = 200) in order not to clutter the 
figure. The frequency of the envelope is e and the dis¬ 
tance between two consecutive valleys is half the period, 
i.e., the distance is l/2|e|. For example, for the systems 
depicted in Figjl] b = V30 = 5.477226, and therefore 
e = —0.022774, resulting in L n ~ 50 for L = 900 and 
L n ~ 54.55 for L = 1000. While for the L = 200 case 
shown in Fig. ^L n ~ 9.1. Indeed, the 9 valley states are 
clearly seen, as well as the fact that the valley positions 
ate not exactly commensurate. 

Thus, to first approximation, there is a superlattice 
of valleys at points j n , with n = l,... ,L n , each with 
a zero-energy state, | n), centered around j = j n . 
These states can be written as Gaussians falling off 
at a length scale of £ = y/f/7rA|e|. The central band 
eigenfunctions are composed of the hybridization of 
these localized state. For periodic boundary con¬ 
ditions, and when L n = \L n \, the eigenstates of 

the central band are plane waves composed of the 
valley Gaussian, | k) = Ln 1 ^ 2 Yln=i S n e lkn \n), where 
S n = \/2cos (nir/2 — 7t/4) = ..., 1,1, —1, —1,1,1 . 
The spectrum Fl central (fc) = —2fcosfc [2(|, 

where k = 2Trm/L n (m = 0, ±1,±2,...) and 

t ft exp(—£ 2 /(4£ 2 e) 2 )(2f exp(—1/4£ 2 ) sinh[(4£ 2 |e|) -1 ] — 
Aexp(—7r 2 e 2 £ 2 )). Thus, the central band spec¬ 
trum is expected to show a degeneracy since 
£central( fc ) = ^central A closer look at this 

issue reveals that if the system is not exactly periodic, 
the degeneracy will be broken by the non perfect 
periodicity. One may think of the effect of the non 
perfect periodicity as an impurity at the region of the 
non-periodicity (i.e., around n = 0). For low-lying 
states in the central band, the impurity acts as a 
hard-wall, leading to low-lying states of the form: 
\k) = \J2/L sin(fcn), with k = 7r m/L n (m = 0,1, 2,...), 
and eigenvalues £' central (fc) = —2tcosfc. Thus, as can be 
seen in the inset of Fig. [TJ for the low-lying states in 
the central band the degeneracy in the eigen-values is 
lost, both for the almost periodic case (L = 900), as for 
the non-periodic one (L = 1000). The low-lying wave 
functions are depicted in the upper panel of Fig. [3] For 


comparison the wave functions of a clean ring of the 
same length, with a single impurity at n = 0 (weak for 
L = 900, strong for L = 1000), is drawn. For the clean 
system the ground state wave function corresponds to 
a half-sine, while the first excited state to a sine. This 
behavior is not very sensitive to the impurity strength. 
A similar situation can be seen for the Harper model, 
where the half-sine and sine envelopes are composed 
of the Gaussian superlattice states at points j n . As 
expected for the low-lying states there is no essential 
difference between the almost periodic case and the 
non-periodic one. 

This changes when higher levels are considered. For 
higher energies, one expects that the effect of the im¬ 
purity will be weaker. Indeed, as can be seen for the 
clean system depicted in the lower panel of Fig. Ql if the 
strength of the impurity is weak (the L = 900 case), the 
wave function is homogeneous. Only for the stronger im¬ 
purity (the L = 1000 case), the wave function shows a 
signature of the impurity. A similar behavior is seen for 
the central band of the Harper model. For the almost 
commensurate case the states in the middle of the cen¬ 
tral band show only a slight degeneracy breaking (see 
the inset of Fig. [I] where the energies for L = 900 
clearly appear in pairs). Nevertheless, when the non- 
commensurability is stronger, as for L = 1000, the levels 
are non-degenerate even in the middle, thus correspond¬ 
ing to a hard wall boundary condition even for higher 
energies. This distinction can be seen also for the wave 
functions, where for the almost commensurate Harper 
model the wave function is almost homogeneous, while 
for the non-commensurate case the wave function is non- 
homogeneous (see Fig. SI). 

Let us now consider the effects of connnensurabilty on 
the persistent current following through a Harper ring. 
The PC of the j- th level is defined as /y(^>) = d£j/dcj)]^, 
while the total persistent current up to a given level n is 
/„((/>) = l g(^)- The total persistent current for any 
value of n for length L = 900 and L = 1000 are presented 
in Fig. [5] with the data pertaining to the energy of the 
n-th level. Gaps between the bands are clearly seen by 
the jumps in the energy between adjacent levels (as in 
Fig. □). First one notes that as expected, the total per¬ 
sistent current reaches a maximum at the middle of the 
bands, while it is zero at the edge of the band. For the 
bands far from middle, there is no essential difference in 
I n between the almost commensurate (L = 900) length 
and the incommensurate length (L = 1000). 

Since we are considering here the metallic regime of 
the Harper model (A = 1), a metallic, behavior of the PC 
is expected. The most obvious difference between the PC 
of a metallic and a localized system is the amplitude of 
the PC. As can be seen from the energy dependence on 
the threading flux presented in Fig. [5j there is a huge dif¬ 
ference in the dependence of the states between almost 
commensurate length ( L = 900 and L = 483) and in- 




FIG. 1: (Color online) The energy spectrums of two different 
length, one with L = 900 (black symbols) and the other with 
L = 1000 (red symbols), for both length b = -\/30 and A = 1. 
The superlattie length correspond to L n ~ 50 for L = 900 
and L n ~ 54.55 for L = 1000. The gross features of the 
spectrum (except for the edge states appearing in the gaps) 
do not essentially change between the commensurate and in¬ 
commensurate length. Inset: a zoom into the central bad for 
both length. Please note the change in the scale of the y axis. 
Here there is a clear difference between the commensurate 



FIG. 2: (Color online) The on-site potential of the Harper 
model for L = 200, with A = 1 and b = \/30. The envelope 
corresponds to cos(27rm). The distance between the valleys is 
l/2e. The number of valleys correspond to L n = 2|e|L which 
for the system depicted in this figure corresponds to L n ~ 9.1, 
which is a slightly incommensurate case. 


commensurate length (L = 1000 and L = 405). In this 
figure the flux dependence of total PC J i / 2 (<)>) is plotted. 
For the commensurate case a high amplitude of the PC is 
seen, which is in line with our expectations from a metal¬ 
lic system. Moreover, a unique sawtooth dependence as 
function of the flux is observed, similar to the typical 
behavior of the total PC in clean and weakly disordered 



FIG. 3: (Color online) The ground and first excited state wave 
functions for L = 900 (upper panel) and L = 1000 (lower 
panel). For the Harper model the ground state wave function 
squared corresponds to the black continuous line, while the 
dashed curve corresponds to the first excited state. For the 
clean ring with an impurity the same notation is used with red 
curves. For the Harper model the Gaussian localized states 
at the superlattice points are seen. It is also evident that the 
Harper states have a similar envelope to the clean state with 
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FIG. 4: (Color online) The L/2+1 excited state wave function 
for the Harper (black curve) and clean system with an impu¬ 
rity (red symbols), for L = 900 (upper panel) and L = 1000 
(lower panel). For L = 900 which is almost commensurate 
(and the corresponding clean system with weak impurity) the 
wave function is more or less homogeneous. For L = 1000 
which is not commensurate (and the corresponding clean sys¬ 
tem with a strong impurity) the wave function shows a spatial 
structure. In both cases the clean system also shows a fast 
(—1)” modulation since it is at half filling. 

systems [23| . On the other hand, for the incommensurate 
case Il/ 2 { 4 > ) is almost flat and has a sine flux dependence, 
expected in the localized regime. 

An additional indicative feature for determining 
whether the PC corresponds to a metallic or a local- 
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FIG. 5: (Color online) The total persistent current for b=\/30 
(e = —0.22774) and A = 1 for different sample length. Upper 
panel: The total PC at a given flux I n (cp = n/2) as function 
of the number of states n for the L = 900 and L = 1000 
length is represented by the black lines. The red curve cor¬ 
responds (with an arbitrary scale) to the energy of the n-th 
level. Jumps in the red curve indicates a gap between the 
states. Close to the gaps (i.e., at the band edge) the total 
PC goes to zero. The main difference between the almost 
commensurate (L = 900) and non commensurate (L = 1000) 
length is the amplitude of the total PC in the central band 
(the states around L/ 2). Lower panel: The total PC at the 
center of the central band as function of the threading flux 
( t>. For the almost commensurate superlattice cases (L = 900, 
and L = 483) the PC is large and it has a saw tooth like shape. 
For the incommensurate cases ( L = 1000, and L = 405) the 
PC is small and the shape is sine like. 


ized behavior, is the number of harmonics composing 
the PC. Generally, the PC may be written as a har¬ 
monic expansion of the form i = t(/) sin(2fmp/ipo) 

and 7 = Hf) sin(2/7r</>/<^ 0 )- For the localized regime 
the PC is described by the first harmonic, i.e., l ~ 
i(l) sin(27r<^/<^o) and 7 ~ 7(1) sin(27r<p/</? 0 ) while 
for the metallic regime, the weight of higher harmonics 


falls off rather slowly [24]. Thus, both the amplitude 
and dependence on the flux of the almost commensu¬ 
rate and non-commensurate length demonstrate a com¬ 
pletely different behavior at the middle of the central 
band. This is also evident from the Fourier tranform of 
the PC (Fig. ©, where for the non-commensurate cases 
only the first harmonic is non-zero, while for almost com¬ 
mensurate cases the contribution of higher harmonics is 
still significant. 

Therefore, if one uses the total PC at half-filling, or 
the persistent current of a particular state in the middle 
of the central band to decide whether the Harper model 
is metallic or localized, one would get different results 
as function of the length of the sample. The behavior is 
metallic for the almost commensurate length and local¬ 
ized like for the non commensurate length. Of course the 
system remains essentially metallic in both cases. Nev¬ 
ertheless, as we have seen from a direct inspection of the 
wave functions in the central band (Figs. 131II) the non 
commensurate segment has an influence similar to an im¬ 
purity embedded in the ring. Since this impurity acts 
like a barrier, we expect that the PC amplitude will be 
suppressed proportionally to the overlap between two su¬ 
perlattice states, i.e., proportionally to exp(C/£ 2 ) (where 
C is a constant). In Fig]7]we plot the first harmonic of 
the PC of the middle state i L i 2 (1) is plotted as function 
of £ 2 = f/7rA|e|. Here we keep the strength of the po¬ 
tential A and system size L constant, while changing |e|. 
It can be seen that ll /2 is suppressed exponentially as 
function of £ 2 . This fits well with the description of the 
non-commensurate central band as a system with an im¬ 
purity. It also confirms that PC can not be used in these 
systems to determine whether the system is metallic or 
localized. 

In conclusion, we have considered the PC in the stag¬ 
gered Harper model, mainly close to the half-filling. Con¬ 
trary to naive expectations the PC for most systems do 
not show typical metallic behavior, although the system 
is in the metallic regime. This stems from the unique 
character of the central band states, which are a hy¬ 
bridization of Gaussian states localized in super lattice 
points. If the distance between the superlattice sites are 
(almost) commensurate with the system size, then the 
PC has metallic properties. On the other hand, when 
the superlattice sites are incommensurate the system has 
an effective impurity at the incommensurate section and 
the PC shows the signature of an insulator. Since for 
most combinations of the Harper frequency b (close to 
the staggered condition Z+l/2) and length L are incom¬ 
mensurate, a typical Harper model will actually exhibit 
PC corresponding to an insulator. 

Thus, studying the PC at half-filling for the finite al¬ 
most staggered Harper model is not indicative of the 
phase of the system. Moreover, also other stout method 
for identifying the metal-insulator transition, such as the 
level spacing statistics turn, will not show the expected 

















5 



f 

FIG. 6: (Color online) Fourier transform of the total PC at 
the middle of the spectrum. For the incommensurate cases 
(L = 1000, and L = 405) the Fourier transform has only 
one significant component which is relatively small. For the 
almost commensurate superlattice cases (L = 900, and L = 
483) the first component is the biggest, but the contribution 
from the higher components goes down only gradually. 

transition in the level statistics from the Poisson distribu¬ 
tion for the localized regime to the Wigner distribution 
in the metallic regime. As is evident from the inset in 
Fig. [TJ the distribution will change even in the metallic 
regime, between a picket fence distribution to a bi-model 
one. Therefore, when studying the Harper model in or¬ 
der to investigate the metal-insulator transition one must 
be careful to choose the appropriate measures. 
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